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ABSTRACT 

We have obtained new spectropolarimetric observations of the planet-hosting star 
t Bootis, using the ESPaDOnS and NARVAL spectropolarimeters at the Canada- 
France-Hawaii Telescope (CFHT) and Telescope Bernard-Lyot (TBL). 

With this data set, we are able to confirm the presence of a magnetic field at the 
surface of r Boo and map its large-scale structure over the whole star. The large-scale 
magnetic field is found to be fairly complex, with a strength of up to 10 G; it features 
a dominant poloidal field and a small toroidal component, the poloidal component 
being significantly more complex than a dipole. The overall polarity of the magnetic 
field has reversed with respect to our previous observation (obtained a year before), 
strongly suggesting that r Boo is undergoing magnetic cycles similar to those of the 
Sun. This is the first time that a global magnetic polarity switch is observed in a star 
other than the Sun; given the unfrequent occurrence of such events in the Sun, we 
speculate that the magnetic cycle period of r Boo is much shorter than that of the 
Sun. 

Our new data also allow us to confirm the presence of differential rotation, both 
from the shape of the line profiles and from the latitudinal shearing that the magnetic 
structure is undergoing. The differential rotation surface shear that r Boo experiences 
is found to be 6 to 10 times larger than that of the Sun, in good agreement with recent 
claims that differential rotation is strongest in stars with shallow convective zones. We 
propose that the short magnetic cycle period is due to the strong level of differential 
rotation. 

With a rotation period of 3.0 and 3.9 d at the equator and pole respectively, r Boo 
appears as the first planet-hosting star whose rotation (at intermediate latitudes) is 
synchronised with the orbital motion of its giant planet (period 3.3 d). Assuming that 
this synchronisation is not coincidental, it suggests that the tidal effects induced by 
the giant planet can be strong enough to force the thin convective enveloppe (though 
not the whole star) into corotation. 

We also detect time dependent activity fluctuations on r Boo, but cannot unam- 
biguously determine whether they are intrinsic to the star or induced by the planet; 
more observations of similar type are needed to determine the role of the close-in giant 
planet orbiting r Boo on both the activity enhancements and the magnetic cycle of 
the host star. 

Key words: stars: magnetic fields - stars: planetary systems - stars: activity - stars: 
individual: r Boo - techniques: spectropolarimetry 



f Based on observations obtained with ESPaDOnS at the Canada-France-Hawaii Telescope (CFHT) and with NARVAL at 
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1 INTRODUCTION 

Magnetic fields of stars hosting close-in giant planets have 
recently started to trigger a lot of interest from the astro- 
physical community. Magnetic fields on the host star are 
indeed likely to play a direct role in the survival of close- 
in giant planets. By evacuating the central regions of their 
protoplanetary discs (w ithi n 0.1 A U from the star typically, 
iRomanova fc Lovelace! e.g., 120061 ) . magnetic fields of young 
protostars are generating ideal conditions for giant planets 
to stop their inward migration whenever they enter the cen- 
tral gap of the disc, leaving them at a location where many 
close-in giant planets (making up 20% of all known exoplan- 
ets to date) are actually observed. 

Magnetic fields of stars hosting close-in giant planets are 
also expected to be a key point in the way such planets in- 
teract with their host stars, and numerous studies have been 
carried out recently to study the details of this interaction. 
Observations suggest that stars hosting close-in giant plan- 
ets exhibit enhanced activity correlating with the orbital 
phase of the planet rather than with the rotation phase of 
the star, arguing that this incr eased activity is induced b y 
the presence of the planet (e.g.. lShkolnik et al.|[20o3 . l2005h ; 
recent observations even suggest that these star-planet in- 
teractions could be cyclic in nature and oscillate between 
'on' and 'off ' states (jShkolnik et al. 2008). The oretical stud- 
ies (e.g. JCuntz et al.ll2000l : iMcIvor et al.ll2006l ) propose that 
this interaction could either result from tidal effects (e.g., en- 
hancing turbulence and hence local dynamo action and ac- 
tivity within the planet-induced tidal bulge) or from magne- 
tospheric interaction (e.g., inducing reconnection events as 
the planet travels through the large magnetospheric loops 
of the host star). 

For such studies, r Boo (HR 5185, HD 120136, F7V) 
is an interesting candidat e. Its massive plan et is orbiting at 
0.049 AU in 3.31 d (e.g., lLeigh et al.ll2003h . Clear Zeeman 
signatures have recently been detected on r Boo with the 
ESPaDOnS spectropolarimeter on the 3.6m Canada- France- 
Hawaii Telescope (CFHT) , demonstrating that a large-scale 
field of a few G is present at the surface (|Catala et al.l 
120071 ). Despite its relatively short rotation p eriod (of or- 
der 3 d. iHenrv et~ai]|2000l ; ICatala et al.l 120071 ). the activity 
of t Boo is only moderate (as usual for mid F stars) and 
shows no strong modulation with either orbital or rotation 
phase (jShkolnik et al.ll2005l .l 20081) ; the near synchronisation 
between the star's rotation and the planet orbit may explain 
this lack of planet-induced activity, r Boo nonetheless ap- 
pears as a good laboratory for studying and modelling the 
magnetic fields of stars with close-in giant planets, and in 
particular those of F stars with shallow convective zones and 
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low intrinsic activity on which the (presumably very small) 
planet- induced activity enhancements are easier to detect. 

We present in the paper a detailed modelling of the 
magnetic field and differential rotation at the surface of 
r Boo, as a follow -up study from the initial discovery of 
ICatala et alj (|2007l ). Sec. [5] presents the observations, Sec. [3] 
details the magnetic and differential rotation modelling 
while Sec.|4]provides a discussion and future prospects about 
this work. 



2 OBSERVATIONS 

For these renewed observations, we used again the ES- 
PaDOnS spectropolarimeter (Donati et al., in preparation) 
on CFHT; we also collected a few additional polarised 
spectra with NARVAL on the 2m Telescope Bernard Lyot 
(TBL). Both instruments yield a spectral resolution of about 
65,000. Each spectrum consists of a sequence of 4 individ- 
ual subexposures taken in different configurations of the po- 
larimeter retarders, in order to perform a precise and achro- 
matic circular polarisatio n analysis and supr ess all spurious 
signatures at first order l|Donati et al.l [19971 . Donati et al., 
in preparation). 

Data were reduced with the dedicated automatic re- 
duction tool Libre- ESpRIT installed at CFHT and TBL 
(IDonati et al. 1997, Donati et al., in preparation), and 
changed into sets of Stokes I and V spectra. All spectra 
are automatically corrected from spectral shifts resulting 
from instrumental effects (eg mechanical flexures, temper- 
ature or pressure variations) using telluric lines as a refer- 
ence. Though not perfect, this procedure allows spectra to 
be secured with a radial velocity (RV) precision of abo ut 
15-20 ms _1 (e.g., [Catala et "al. I l2007l : iMoutou et al.|[2007l ). 

A total of 32 spectra were collected in 2007 June 
and July in variable weather conditions, mostly with ES- 
PaDOnS / CFHT. The complete log is given in Table [TJ All 
data are pha s ed wi th the same orbital ephemeris as that of 
ICatala et all (|2007m : 



To = HJD 2, 453, 450.984 + 3.31245S 



(1) 



with phase 0.0 denoting the first conjunction (i.e., with the 
planet farthest from the observer). Circularly polarised spec- 
tra of stars with stable magnetic topologies (e.g., chem- 
ically peculiar sta rs and hotter equivalents, like r Sco, 
IDonati et ai1 l2006) were observed during both runs to check 
that the instrument behaves properly and yields nominal 

field strengths and polarities. 

Least-Squares Deconvolution (LSD, IDonati et al.lll997h 
was applied to all spectra to retrieve average unpolarised and 
circularly polarised profiles of photospheric spectral lines. 
The line pattern used for this process is derived from a 
Kurucz model atmosphere with solar abundances, and ef- 
fective temperature and logarithmic gravity set to 6250 K 
and 4.0 respectively; this line pattern includes most mod- 
erate to strong lines present in the optical domain (those 
featuring central depths larger than 40% of the local contin- 
uum, before any macroturbulent or rotational broadening, 
about 4,000 throughout the whole spectral range) but ex- 
cludes the very strongest, broadest features, such as Balmer 
lines, whose Zeeman signature is strongly smeared out com- 
pared to those of narrow lines. The typical multiplex gain 
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Table 1. Journal of observations. Columns 1-6 sequentially list the UT date, the instrument used, the heliocentric Julian date and 
UT time (both at mid-exposure), the complete exposure time and the peak signal to noise ratio (per 2.6 kms" 1 velocity bin) of each 
observation. Column 7 lists the rms noise level (relative to the unpolarized continuum level I c and per 1.8 kms -1 velocity bin) in the 
circular polarization profile produced by Least-Squares Deconvolution (LSD), while column 8 and 9 list the orbital cycle (using the 
ephemeris given by Eq. [TJ and the radial velocity (RV) associated with each exposure. 



Date Instrument HJD UT i oxp S/N <?lsd Cycle ^ rad 

(2007) (2,453,000+) (h:m:s) (s) (10 -4 / c ) (245+) (kms" 1 ) 



Jun 


12 


NARVAL/TBL 


1264. 


.45160 


22: 


:46: 


:05 


4 


X 


600 


1,000 


0. 


37 


0. 


.5788 


-16 


.506 


Jun 


19 


NARVAL/TBL 


1271. 


.42056 


22: 


:02 


:04 


4 


X 


600 


910 


0. 


12 


2. 


.6827 


-16 


.708 


Jun 


26 


ESPaDOnS / CFHT 


1277. 


.74494 


05 


:49: 


:50 


4 


X 


300 


1,550 


0. 


23 


4. 


.5920 


-16 


.589 


Jun 


26 


ESPaDOnS / CFHT 


1277. 


.81018 


07 


23: 


:47 


4 


X 


300 


1,730 


0. 


18 


4. 


.6117 


-16 


.645 


Jun 


27 


ESPaDOnS / CFHT 


1278. 


.74598 


05 


:51 


:27 


4 


X 


200 


1,700 


0. 


21 


4. 


.8942 


-16 


.644 


Jun 


27 


ESPaDOnS / CFHT 


1278. 


.75739 


06: 


:07 


:53 


4 


X 


200 


1,670 


0. 


21 


4. 


.8976 


-16 


.641 


Jun 


27 


ESPaDOnS / CFHT 


1278. 


.76861 


06: 


:21 


:02 


1 


X 


200 


1,710 


0. 


21 


4. 


.9010 


-16 


.625 


Jun 


28 


ESPaDOnS / CFHT 


1279. 


.73666 


05: 


:38: 


:08 


4 


X 


200 


1,620 


0. 


22 


5. 


.1933 


-15 


.927 


Jun 


28 


ESPaDOnS / CFHT 


1279. 


.74787 


05 


:51 


:17 


4 


X 


200 


1,660 


0. 


22 


5. 


.1966 


-15 


.915 


Jun 


28 


ESPaDOnS / CFHT 


1279. 


.75913 


06: 


:10 


:30 


4 


X 


200 


1,730 


0. 


21 


5. 


.2000 


-15 


.912 


Jun 


30 


ESPaDOnS / CFHT 


1281. 


.83298 


07: 


:57. 


:04 


4 


X 


200 


1,470 


0. 


25 


5. 


.8261 


-16 


.746 


Jun 


30 


ESPaDOnS / CFHT 


1281. 


.84420 


08. 


:13 


:13 


4 


X 


200 


1,290 


0. 


28 


5. 


.8295 


-16 


.739 


Jul 


01 


ESPaDOnS / CFHT 


1282. 


.73612 


05: 


:37 


:41 


1 


X 


200 


1,640 


0. 


21 


6. 


.0988 


-16 


.070 


Jul 


01 


ESPaDOnS /CFHT 


1282. 


.74732 


05 


:53 


:49 


1 


X 


200 


1,640 


0. 


21 


6. 


.1022 


-16 


.062 


Jul 


01 


ESPaDOnS / CFHT 


1282. 


.75863 


06 


:10 


:06 


4 


X 


200 


1,720 


0. 


21 


6. 


.1056 


-16 


.058 


Jul 


02 


ESPaDOnS / CFHT 


1283. 


.73836 


05: 


:41: 


:01 


4 


X 


200 


740 


0. 


11 


6. 


.4013 


-16 


.093 


Jul 


02 


ESPaDOnS / CFHT 


1283. 


.74973 


05 


:57. 


:24 


4 


X 


200 


580 


0. 


53 


6. 


.4048 


-16 


.090 


Jul 


02 


ESPaDOnS / CFHT 


1283. 


.76769 


06: 


:23 


:16 


4 


X 


200 


1,240 


0. 


31 


6. 


.4102 


-16 


.103 


Jul 


02 


ESPaDOnS / CFHT 


1283. 


.78209 


06: 


:43: 


:60 


1 


X 


200 


1,360 


0. 


27 


6. 


.4145 


-16 


.110 


Jul 


02 


ESPaDOnS / CFHT 


1283. 


.91960 


10 


:02 


:02 


4 


X 


300 


1,370 


0. 


29 


6. 


.4561 


-16 


.193 


Jul 


03 


ESPaDOnS / CFHT 


1284. 


.74158 


05: 


:45: 


:46 


4 


X 


300 


1,370 


0. 


24 


6. 


.7042 


-16 


.816 


Jul 


0.3 


ESPaDOnS / CFHT 


1284. 


.75763 


06: 


:08. 


:53 


4 


X 


300 


1,830 


0. 


19 


6. 


.7090 


-16 


.818 


Jul 


0.3 


ESPaDOnS / CFHT 


1284. 


.77363 


06: 


31: 


:56 


4 


X 


300 


1,750 


0. 


21 


6. 


.7139 


-16 


.808 


Jul 


03 


ESPaDOnS / CFHT 


1284. 


.85319 


08. 


:26 


:30 


4 


X 


300 


1,700 


0. 


23 


6. 


.7379 


-16 


.817 


Jul 


04 


ESPaDOnS / CFHT 


1285. 


.79220 


06 


:58: 


:47 


4 


X 


300 


1,090 


0. 


33 


7. 


.0214 


-16 


.275 


Jul 


04 


ESPaDOnS / CFHT 


1285. 


.80803 


07: 


:21 


:35 


4 


X 


300 


1,470 


0. 


26 


7. 


.0262 


-16 


.248 


Jul 


04 


ESPaDOnS / CFHT 


1285. 


.86573 


08: 


:44: 


:40 


4 


X 


300 


1,240 


0. 


32 


7. 


.0436 


-16 


.191 


Jul 


04 


ESPaDOnS / CFHT 


1285. 


.88171 


09 


:07 


:41 


4 


X 


300 


990 


0. 


41 


7. 


.0484 


-16 


.173 


Jul 


05 


ESPaDOnS / CFHT 


1286. 


.74025 


05 


:44: 


:05 


4 


X 


200 


330 


1. 


18 


7. 


.3076 


-15 


.943 


Jul 


05 


ESPaDOnS / CFHT 


1286. 


.75144 


06: 


:00 


:12 


1 


X 


200 


250 


1. 


.57 


7 


.3110 


-15 


.933 


Jul 


05 


ESPaDOnS / CFHT 


1286. 


.76265 


06: 


:16 


:20 


4 


X 


200 


310 


1. 


28 


7. 


.3143 


-15 


.930 


Jul 


05 


ESPaDOnS / CFHT 


1286. 


.78415 


06: 


:47: 


:18 


4 


X 


600 


400 


0. 


98 


7 


.3208 


-15 


.941 



in S/N for polarisation profiles is about 25 to 30, implying 
noise levels in LSD polarisation profiles as low as 20 ppm 
(i.e., 2 X 10 - in units of the unpolarised continuum) in good 
weather conditions. Zeeman signatures, whenever detected, 
have a typical amplitude of about 100 ppm (see below). 

Radial velocities are obtained from Gaussian fits to each 
LSD unpolarised profile of r Boo. We find that our measure- 
ments are in good agreement with previous measurements 
(see Fig. [1} . Residuals (rms) with respect to the predicted 
RV curve are equal to 20 ms -1 , i. e., similar to the u ncer- 
tainties of 15-20 ms " 1 reported bv lCatala et alj (|2007f ) and 
iMoutou et~aH (|2007T ). All Stokes / and V profiles used in the 
following were corrected for the orbital motion, i.e., trans- 
lated into the velocity rest frame of r Boo. 

As in ICatala et al. 

I (|2007l ). we assume that r Boo ro- 
tates about an axis inclined at an angle i = 40° with re- 
spect to the line-of-sight. The unpolarised spectral lines are 
significantly broadened by rotation and suggest that the 
rotation of the star is more or less synchronised with the 
planet orbital motion, i.e., corresponds to a rotation rate 
of 1.9 radd -1 . By averaging all Stokes / LSD profiles into 
a single mean line and computing its Fourier transform, 



one can estimate th e amount of differential r otation at the 
surface of the star (Rei ners fc Schmittll2002l ). This experi- 
ment, first c a rried out by iReiners! (|2006l ) and repeated by 
ICatala et all i|2007l ), demonstrates that r Boo is indeed dif- 
ferentially rotating, with a relative differential rotation (i.e., 
an angular rotation shear relative to the mean angular ro- 
tation rate) of 18-20% and an angular rotation sheaiQ of 
0.35 ± 0.10 radd -1 ; the corresponding rotation periods at 
the equator and the pole are equal to about 3.0 d and 3.7 d 
respectively, bracketing the planet orbital period of 3.3 d. 
The Fourier transform of the average LSD Stokes / profile 
derived from our ne w data (not shown here) is very simi- 
lar to that shown in lCatala et al.l l|2007l ) (their Fig. 1) and 
confirms their analysis. 



1 The angular rotation shear derived by I Reiner 1 (I2006T) . equal 
to 0.31 ± 0.13 radd -1 , is actually scaled down by V sin i; using 
i = 40°, it translates into a true angular rotation shear of 0.38 ± 
0.16 radd -1 in good agreement with our own estimate. 
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1280 1285 0.2 0.4 0.6 0.8 

Heliocentric Julian Date (2,453,000+) Orbital phase 



Figure 1. Radial velocities of r Boo derived from our ESPaDOnS spectra (green circles) as a function of Heliocentric Julian Date (left) 
and orbital phase (right, using the ephemeris of eq. [TJ. The ra dial velocity model plotted here (red full line) corresponds to a circular 
orbit with an (assumed) velocity semi-amplitude of 467 ms~ _1 l|Butler et al.lll997h and a (fitted) systemic velocity of —16.348 kms . 
Individual error bars on data points (15-20 ms _1 ) are about as large as the symbol size. 



3 MAGNETIC MODELLING 
3.1 Model description 

To model the m agnetic topology of r Boo, we use the new 
imaging code of IDonati et al.l (|2006l ) where the field topol- 
ogy is described through spherical-harmonics expansions. 
We use the principles of maximum entropy image recon- 
struction to retrieve the simplest magnetic image compatible 
with the series of rotationally modulated Zeeman signatures. 
More specifically, the field is divided into its radial-poloidal, 
non-radial-poloidal and toroidal components, each of them 
described as a spherical-harmonics expansion; given the 
different rotational modulation of Zeeman signatures that 
poloidal and toroidal fields generate, our imaging code ap- 
pears particularly useful and efficient at producing dynamo- 
relevant diagnostics about the large-scale magnetic topolo- 
gies at the surface of late-type stars. One of the latest a p- 
plication of this code can be found in iMorin et all (|2007l ). 

The reconstruction process is iterative and proceeds by 
comparing at each step the synthetic profiles corresponding 
to the current image with those of the observed data set. 
To compute the synthetic Zeeman signatures, we divide the 
surface of the star into small grid cells (typically a few thou- 
sands) , work out the specific contribution of each grid cell to 
the Stokes V profiles (given the magnetic field strength and 
orientation within each grid cell, as well as the cell radial ve- 
locity, location and projected area) and finally sum up con- 
tributions of all cells. Since the problem is partly ill-posed, 
we stabilise the inversion process by using an entropy crite- 
rion (applied to the spherical harmonics coefficients) aimed 
at selecting the image with mini mum information a mong all 
those compatible with the data l|Morin et al.ll2007l ). 

The model we use to describe the local Stokes I and V 
profiles is quite simple. While the local unpolarised profile 
is given by a Gaussian, the local circular polarisation pro- 
file is computed assuming the weak field approximation, i.e., 
that V is proportional to dl /d\ and to the local line-of-sigh t 
component of the magnetic field (e.g.. IDonati et ai1ll997h . 
Assuming a line-of-sight projected equatorial rotation veloc- 
ity «sini of 15.9 kms -1 and a relative differential rotation 



of 18%, we obtain synthetic Stokes / profiles whose first 2 
zeros in the Fourier transform match those of the observed 
data; setting the full width at half maximum of the local 
Stokes / profile to 11 kms -1 produces a very good fit to the 
average Stokes I LSD profiles. This simple line model was 
used quite extensively and has proved to be efficient at cor- 
rectly reproducing obser ved sets of rotation a lly modulated 
Zeem an signatures (e.g., ICatala et al. 1 120071 ; iMoutou et al.l 
120071) . 

To incorporate differential ro t ation into t he modelling, 
we pr oceed as in IDonati et ail (|2003l ) and |m orin et al.l 
|2007l ). i.e., by assuming that the rotation rate at the sur- 
face of the star is varying with latitude as fi eq — dQ sin 2 6 
where f2 Gq is the rotation rate at the equator dfl the dif- 
ference in rotation rate between the equator and the pole. 
When computing the synthetic profiles, we use this law to 
work out the longitude shift of each cell at each observing 
epoch with respect to its location at the median observing 
epoch (at which the field is reconstructed, i.e., orbital cycle 
6 + 245 = 251 in the present case and in the ephemeris of 
eq.[T]or HJD 2,454,282.41) so that we can correctly evaluate 
the true spectral contributions of all cells at all epochs. 

For each pair of f2 eq and d£l values within a range of 
acceptable values, we then derive, from the complete data 
set, the corresponding magnetic topology (at a given in- 
formation content) and the associated reduced chisquare 
level x? a t which modelled spectra fit observations. By fit- 
ting a paraboloid t o the Xt surface derived in this process 
l|Donati et al.|[200 3). we can easily infer the magnetic topol- 
ogy that yields the best fit to the data along with the cor- 
responding differential rotation parameters and error bars. 
This process has proved fairly reliable to est imate surface 
differ ential rotation on magnetic stars (e.g., IDonati et al.l 
I2003D . 

3.2 Results 

We applied this model to our complete set of Zeeman sig- 
natures of r Boo. Data obtained on the last night (Jul 05), 
i.e., in bad weather conditions, were finally excluded from 
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2 2.05 2.1 2.15 2.2 



eq (rad/d) 

Figure 2. Variation of Xt as a function of f2 eq and dfi, derived 
from the modelling of our Stokes V data set on r Boo. The outer 
colour contour corresponds to a 1% increase in x?> an d trace 
a 2.7a interval for each parameter taken separately, or a 1.7a 
interval for both parameters as a pair. 



the fit as they provide very little information (given their 
low quality with respect to the bulk of the data set, see Ta- 
ble [TJ. The rotational broadening of spectral lines provides 
significant spatial resolution at the surface of the star, i.e., 
from 9 to 22 resolution elements around the equator depend- 
ing on whether we use the intrinsic line profile (11 kms -1 ) 
or the instrumental profile (4.5 kms -1 ) to define the res- 
olution clement. The spherical harmonics expansions used 
to describe the magnetic field were first computed up to or- 
ders I = 15; in practice, little improvement is obtained when 
adding terms with orders higher than 8 in the expansion. All 
results presented below were thus derived with expansions 
limited to £ = 8. 

The first result from this modelling is that the Stokes 
V data provide a completely independent confirmation of 
both the existence and magnitude of differential rotation 
at the surface of the star, as demonstrated by the well de- 
fined x? paraboloid we obtain (see Fig. [2}. The differen- 
tial rotation parameters producing an optimal fit to the 
data (at a given information content in the magnetic map) 
are respectively equal to Q cq = 2.10 ± 0.04 radd -1 and 
dfi = 0.50 ± 0.12 radd -1 and imply a relative differential 
rotation of 24%. The corresponding rotation periods at the 
equator and pole are equal to about 3 d and 3.9 d respec- 
tively; the optimal period assuming solid body rotation is 
3.23 d and corresponds to the average recurrence time of 
the detected Zeeman signatures, i.e. to the rotation period at 
the average latit ude (about 35°) of the reconstructed mag- 
netic structures l|Donati et al.l l 2003). The latitudinal surface 
shear we derive is about 8 to 10 times stronger than that of 
the Sun; the corresponding timescale for the equator to lap 
the pole by one complete rotation cycle is only about 12 d. 

The large-scale magnetic map we derive (see Fig. [3]) 
corresponds to a unit Xv fit to the data (see Fig. |4j, the 
initial Xr being about 1.9 for 750 data points. Note that 
the evolution of the magnetic field under the shearing effect 
of differential rotation throughout the period of our obser- 
vations can be traced directly to the data themselves; for 



instance, the Zeeman signature collected at cycle 2.68+245 
is significantly smaller in amplitude than those collected 4 
rotation cycles later (around cycles 6.70+245). 

We find that the field topology includes a small toroidal 
component enclosing 17% of the overall reconstructed mag- 
netic energy; while the Stokes V data can be fitted without 
the toroidal component, the corresponding reconstructed 
map contains significantly more (i.e., +50%) information 
than that of Fig. [3] suggesting that a purely poloidal field 
is less probable and that the reconstructed toroidal compo- 
nent is likely real. The toroidal component is clearly visible 
in Fig.|3]and shows up as ring of positive (i.e., counterclock- 
wise) azimuthal field encircling the star at mid latitude. 

The magnetic field is obviously more complex than a 
dipole, 1=1 modes enclosing only about 30% of the poloidal 
field energy. The quadrupolar and octupolar terms (i.e. 1 = 2 
and I = 3 modes) dominate the field distribution in the 
visible hemisphere and contain 40% of the reconstructed 
poloidal field energy, while the remaining 30% spreads into 
higher order terms. This is directly visible from Fig. [3] the 
radial field map features a main positive pole at high lati- 
tudes surrounded by an incomplete ring of negative field at 
low latitudes, reminiscent of a slightly tilted quadrupole or 
octupole. We also find that the reconstructed poloidal field 
is mostly axisymmetric with respect to the rotation axis; 
modes verifying m < 1/2 are enclosing a dominant fraction 
(i.e. 60%) of the poloidal field energy, while the non axisym- 
metric modes (with m > 1/2) contain no more than 30% of 
the poloidal field energy. This is again fairly obvious from 
the reconstructed map of Fig. [3] 



4 DISCUSSION AND CONCLUSION 

Thanks to this new data set, we achieved a number of sig- 
nificant results relevant to dynamo processes and magnetic 
field generation in cool stars with shallow convective zones; 
We summarise them below and discuss their implications 
for the study of star/planet magnetic interactions in sys- 
tems hosting close-in giant planets. 

First, we obtained 2 completely independent estimates 
of the differential rotation at the surface of r Boo. From the 
temporal distortion of the large-scale magnetic topology, we 
find that the latitudinal angular rotation shear is equal to 
dQ. = 0.50 ± 0.12 radd" 1 , i.e., about 8 to 10 times that of 
the Sun. This is in reasonable agreement with the estimate 
derived from the detailed shape of spectral lines and of their 
Fourier transform, yielding dQ, = 0.35 + 0.10 radd - . It un- 
ambiguously demonstrates that r Boo is experiencing strong 
differential rotation at photospheric level; th is is apparently 
a general trend of early G and F stars fe.g.. iMarsden et al.l 
|2006| ;Ir einersl l2006). Our measurement is also the first direct 
and simultaneous confirmation that both methods employed 
up to now to investigate latitudinal shears on stellar surfaces 
are actually yielding consistent results. The rotation period 
at the equator of r Boo is 3 d, i.e., about 10% shorter than 
the orbital period of the giant plane t. It confirms in par tic- 
ular the estimate first obtained by ICatala et alj (|2007T l; it 
also implies that the giant planet is synchronised with the 
surface of the star at a latitude of about 40° . 

We also derived how the magnetic field is distributed at 
the surface of r Boo. In particular, we find that the mag- 
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Radial magnetic field 
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Azimuthal magnetic field 
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Figure 3. Maximum-entropy reconstructions of the large-scale magnetic topology of r Boo as derived from our 2007 data set (at orbital 
cycle 6+245 or HJD 2,454,282.41). The radial, azimuthal and meridional components of the field are displayed from left to right (with 
magnetic flux values labelled in G). The star is shown in flattened polar projection down to latitudes of —30°, with the equator depicted 
as a bold circle and parallels as dashed circles. Radial ticks around each plot indicate orbital phases of observations. 
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Figure 4. Maximum-entropy fit (thin red line) to the observed Stokes V LSD profiles (thick black line) of r Boo. The orbital cycle of 
each observation (as listed in Table 1) and la error bars are also shown next to each profile. 



netic field is mostly poloidal despite t he vigorous differentia l 
rotation. This is different than what iMarsden et all (|2006l ) 
report for another cool star with a shallow convective zone 
(HD 171488), in which the magnetic field is apparently dis- 
tributed roughly evenly between poloidal and toroidal field 
components; the main difference between both stars is the 



rotation rate, less than half as large for r Boo than for 
HD 171488. The magnetic topology we reconstruct is grossly 
si milar to (though mu ch more accurate than) that derived 
bv lCatala et al.l (|2007l ) from a much sparser data set. There 
is however one major difference between both maps; while 
the radial field is predominantly positive and negative at 
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high and low latitudes respectively in our ima ge, the oppo- 
site h olds in the 2006 image reconstructed bv lCatala et al.l 
(2007). A similar polarity inversion is observed for the two 
other field components between the 2006 and 2007 images. 
Check stars with known magnetic polarities observed dur- 
ing both runs demonstrate that this global polarity switch 
is not due to instrumental or data reduction problems and 
can only be attributed to r Boo itself. 

This is the first time that a global magnetic polarity 
switch is observed in a star other than the Sun. Given that 
such events are rather unfrequent in the Sun (only once every 
11 yr) and have never been detected yet in the 20 or so stars 
(of various spectral types ) that have been observed more 
than once up to now (e.g.. iDonati et aTll2003l ). we speculate 
that the magnetic cycle of r Boo is likely shorter than that of 
the Sun. If the cycle period varies more or less linearly with 
the strength of differential rotation, we expect the period 
of the full magnetic cycle (22 yr in the case of the Sun) 
to be of order 2-3 yr in r Boo; this argues for renewed 
and regular spectropolarimetric observations of r Boo to 
monitor the evolution of the magnetic field throughout a 
complete magnetic cycle. Given that the convective zone of 
r Boo is very shallow and essentially reduces to a thin layer 
similar in nature to the solar tachocline (where the magnetic 
dynamo is expected to operate mostly), our result directly 
demonstrates that interface dynamos are indeed capable of 
producing oscillating magnetic topologies. 

At this stage, there is not much we can say about activ- 
ity putatively induced by the presence of the giant planet. 
There is actually no low-latitude magnetic features either 
facing the planet (at phase 0.5) or on the other side of the 
star (phase 0.0), i.e., in regions at which tidal effects are 
maximum. The main low-latitude magnetic features we de- 
tect are the negative radial and meridional field features at 
phases 0.13, 0.40, 0.60 and 0.93 (see Fig. [3}; however, these 
features are apparently rotating faster than the planet or- 
bital motion (being those from which differential rotation at 
low latitudes is estimated) and can therefore not be inter- 
preted as due to a putative tidal bulge (rotating in phase 
with the orbital motion). 

We find that activity signatures in usual spectral in- 
dexes (Ha, Ca II H & K, and infrared triplet lines) are very 
weak, smaller than 0.5% of the unpolarised continuum; they 
are best visible in Ha thanks to the higher spectrum quality 
around 700 nm (see Fig. [SJ) and reach a peak to peak am- 
plitude of about 0.2 kms -1 (0.44 pm) only. Maximum Ha 
emission (i.e. positive integrals over the residual spectra of 
Fig. [SJ) occurs roughly twice per rotation (at phases 0.1 and 
0.7) and roughly coincide with the main low-latitude radial 
field features seen in the magnetic map; we therefore suspect 
that the two small activity enhancements we detect trace in- 
trinsic activity from the star itself. It is unlikely that this 
activity is induced by the planet itself through tidal fric- 
tion (also expected to produce two activity enhancements 
around phase 0.0 and 0.5); given that the stellar equator is 
rotating faster than the orbital motion, the tidal bulge is 
expected to be slightly ahead of the planet (shifted to neg- 
ative orbital phases), i.e., the opposite of what we actually 
see. The near synchronisation between the planet orbital mo- 
tion and the star equatorial rotation (beat period of about 
32 d) likely implies that the tidal bulge generates very lit- 
tle friction and features only a marginal misalignment with 
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Figure 5. Residual Ha signature as a function of orbital cycle. 
Profiles recorded on the same night were averaged together to 
reduce noise, and a mean Ha profile was subtracted from each 
nightly average. This operation was performed on ESPaDOnS 
spectra only. The stray pixels at +50 kms - are due to a weak 
(and variable) telluric line. 



the star-planet direction. Magnetic reconnection triggered 
by the planet nevertheless remains a potential option for 
explaining the observed activity variations. 

We note that r Boo is apparently the first convincing 
case of a star whose rotation is synchronised wi th the orbita l 
mot ion of its clo s e-in g iant planet. Following IZahnl (|l994l ) 
and iMarcv et all (|l997l ). the synchronisation timescales t s 
(in yr) it takes for the planet to enforce corotation of the 
whole star can be approximated by: 



t s = 4.5 



(2) 



where M/m is the star to planet mass ratio (equal to about 
350 sin i for r Boo) and a/R the planet orbit semi-major 
axis relative to the radius of the star (equal to about 7.2 for 
t Boo). Using i — 40°, we obtain that t s ~ 30 Gyr, i.e., much 
larger than the estimated lifetime of r Boo (about 1 Gyr). 
Since we consider that the apparent spin-orbit synchronisa- 
tion between r Boo and its close-in giant planet is unlikely 



8 J.-F. Donati et al. 



to be coincidental, we suspect the synchronisation timescale 
to be likely overestimated; one possible reason is that syn- 
ch ronisation is not achieved on the whole star (as assumed 
bv lZahnlll994r ) but only on a restricted volume close to the 
surface of the star, e.g., the convective zone proper whose 
mass is estimated to less than 0.1% that of the star. If the 
tidal effects induced by the close-in giant planet on its host 
star are strong enough to trigger synchronisation of at least 
the shallow convective zone, they may also play a significant 
role in the dynamo processes operating in this thin layer. 

More high-quality observations densely sampling the or- 
bital and rotation cycles, carried out over a timescale of 
typically a month, and repeated at least once a year, are re- 
quired to go further along these tracks. Such data sets will 
first allow us to obtain a complete magnetic monitoring of 
the activity cycle of r Boo and give us the opportunity to 
achieve the first such study on a star other than the Sun. 
These data should also enable us to estimate the lifetime 
and recurrence rate of the activity enhancements that we 
detected on r Boo; if these episodes turn out to be long- 
lived and synchronised with the orbital motion (rather than 
with the rotation rate at the stellar equator), they could 
be unambiguously attributed to the giant planet. Finally, 
carrying out similar observations on a sample of stars with 
and without planets to look for statistical differences be- 
tween both subsamples will also be necessary to investigate 
in more details the impact of close-in giant planets on the 
dynamo processes of stars with shallow convective zones. 
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